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ff^ ' Abstract 

o' 

0^ ' The entropy-to-energy bound is examined for a quantum scalar field confined to a 

o ■ 

O . 

^^ , cavity and satisfying Robin condition on the boundary of the cavity. It is found that 

"o 

O^' near certain points in the space of the parameter defining the boundary condition 

rV' the lowest eigenfrequency (while non-zero) becomes arbitrarily small. Estimating, 

^ . following Bekenstein and Schiffer, the ratio S/E by the (^-function, (24(^(4)) ', we 

compute C(4) explicitly and find that it is not bounded near those points that signals 
violation of the bound. We interpret our results as imposing certain constraints on 
the value of the boundary interaction and estimate the forbidden region in the 
parameter space of the boundary conditions. 
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1 Introduction 

Some time ago Bekenstein proposed ^ , that for a quantum system confined to a cavity 
of finite size R the entropy to the energy ratio S/E can not be arbitrarily large so that 
the bound 

S/E < 27rR (1) 

takes place. Originally the bound was deduced by considering a gedanken experiment of 
lowering the system into a black hole and demanding this process to satisfy the generalised 
second law. The bound (|ID comes out as a consistency condition between the black hole 
thermodynamics and ordinary statistical physics. Since the system initially can be placed 
far away from the black hole where the gravitational field is negligible, the bound must 
hold for any system in flat space-time and be provable with no recourse to gravitational 
physics. Thus, the universality of (|l|) was conjectured. 

The black hole way of deriving the bound was criticised in (see, however [Q). 
Nevertheless, the bound (|l]), as it stands, has passed a number of tests ||l|, 0, 0, 0, 
1^ (for a review see [§) so that its universality deserves the further examination (for a 
recent discussion see and [0). In order to make the statement on the entropy bound 
precise one has to define the meaning of S, E and R in (|1]) as well as the conditions under 
which the statistical properties of the system should be considered. In ref. the bound 
(|I]) is regarded as applying only to the field in the cavity and it is proposed to use the 
microcanonical methods. One interprets S as logarithm of ^{E), the number of quantum 
states accessible to the field system with energy up to and including E, and ignores the 
walls of the cavity. 

The bound (|I|) can be obviously b exceeded if there are one-particle states with zero- 
energy (zero- modes) ||rT|. Then, by adding arbitrary number of such states one does 
not change the total energy of the system but makes the entropy S arbitrary large. 



The important observation |1^, however, is that the zero- mode with some occupation 
corresponds to a condensate. The systems with different configurations of the condensate 
should be considered as macroscopically different. Thus, only the excitations with energy 
above the vacuum should be taken into account, i.e. the zero-modes are to be excluded. 



Suppose, the cavity confining the system is circumscribed by a sphere of radius R. 
Then, it was shown in [Q that the microcanonical entropy S{E) = lnQ{E) obeys 

S{E)/E < [24Cp(4)]V^ , (2) 

where (sp{k) is the ^-function 



Cik) = E 



UJ, 



for the sphere, where {ui} is the discrete one-particle energy spectrum with zero-modes 
excluded. Since for the sphere we have ({4) ~ R^ the bound (0) follows from (^ provided 
R~'^({A) is appropriately bounded from above. The later was verified in [^ for various 
types of free fields satisfying Dirichlet or Neumann conditions on the sphere. 

In this paper we make a step further and impose more general, of Robin type, condition 
on the field on the boundary of the cavity. Note, that the boundary condition of this type 
should be always imposed on a quantum field non-minimally coupled to the metric. The 
simplest case is the scalar field described by the action 



2 Jm ^ ' JdM 



where R is Ricci scalar and K is the extrinsic curvature of the boundary DM. The 
boundary term is necessary to add in order to the stress-tensor for the theory be well- 
defined. Variation of this action with respect to (j) gives us not only the equation of motion 
in the interior 

-V^0 + iR(t) = Q 

but also the boundary condition 

Ka^0 + 2eir0),^ = O (3) 

of the Robin type. Note also that the allowing for a more general boundary condition 
is in accord with the general assumption of the Schiffer-Bekenstein paper that all 
interactions of the field are negligible "except for those which confine it and are expressed 
as boundary conditions". The boundary condition (^ encodes in a generic form such 
boundary interaction ||13|. 



The (^-function is simpler object for computation (the computational technique ap- 
propriate to the case under consideration was developed in [jl4| and[0) than the ratio 
S{E)/E. Therefore, in this paper we mainly analyse the (^-function. Considering R~^({A) 
as function of the parameter in the Robin boundary condition we find that there are 
special points in the parameter space near which this function is unbounded from above. 
This can be easily understood. Exactly at those special points the quantum field has a 
zero-mode. When any of these points is approached in the parameter space it signals in 
that the lowest eigenfrequency uJi, while non-zero, becomes arbitrarily small. Since Ui 
makes the dominant contribution to the ^-function (see also 0), one has (^(4) ~ 4r if the 
lowest energy state is non-degenerate. It is evident that C(4) is unbounded in this case. 

The same is also true for the ratio S{E)/E itself. When the parameter in the boundary 
condition approaches one of those special points it is typical that in the spectrum appears a 
large gap between the lowest (non-zero) eigenfrequency ui and the next eigenfrequency uj2, 
UJ2/UJ1 >> 1. When the energy of the system is E = nEi for some integer n and E < uj2 
only the lowest energy level is populated. The number of accessible states (assuming, for 
simplicity, that gi = 1 for the degeneracy of the lowest energy level) is VL{E) = {n + 1) 
and we have S{E)/E = —n^^ln{n+ 1). Since max(?2~-'^ ln(?2 + 1)) = ln2 ~ 0.7 we obtain 
that [|1|] 

max(S(E)/E) ~ — . (4) 

This shows, in particular, that relation (^ should be considered as a good estimate for 
the maximum of S{E)/E rather than just giving an upper bound on S{E)/E. It follows 
that the bound (|I]) holds only if Rui is restricted from below and is violated if Rui can 
be made arbitrarily small. The later does occur for the certain values of the parameter 
in the boundary condition as we show in this paper. 

In the next section we consider in detail the case of (1 + l)-dimensional massless field 
for which the analysis of the spectrum and the computation of (^-function are especially 
simple. The (3 + 1) -dimensional field is analysed in section 3 and the massive field is 
briefly discussed in section 4. The size of "forbidden region" in the parameter space is 
estimated in section 5. Section 6 contains some concluding remarks. 



2 The interval 



We start our analysis with consideration of (1 + l)-diniensional massless scalar field living 
on the interval [0, R] with the boundary conditions 



dx R 







x=0 



— < 

dx 



h 
r' 







(5) 



x=R 



The energy eigenfunction 0(^(x)e*^* satisfies the differential equation 

d 



and takes the form 

4>u;{x) = N sm{ujx + 6) . 

The spectrum is discrete, Un = anR~^, where {a} are (positive) roots of the equation 

tan a ho + hi 
a hohi — a^ 

In general, there also may be bound states (for which uJ^ - 
function (pxlx) = AeR^ + Be^R^. The equation on A is 

tanh A ho + hi 



(6) 



(7) 



-A^i? ^ < 0)) with wave 



A 



A^ — hohi 



As a consequence of the Mittag-Leffler theorem, we have |T5 

Sill z 

{hohi - z^) {ho + hi) cosz = {hohi -ho- hi) W{'\- - 



a^ 



(9) 



This formula helps to evaluate explicitly the sums of inverse powers of the roots. One 
just has to take the logarithm of both sides of eq.(pl), expand in powers of z^ and equate 



the relevant coefficients. The expressions for any ho and hi are given in |^. To simplify 
the things, in what follows, we assume that ho = 0,hi = h. One finds. 



V— -1-1 



a^ 



h 



^^ ^,4 U2 QU R 



a" 



/i2 3h Q 



(10) 



:ii) 



Let us consider the case of positive h first. For ho = 0,hi = h the equation (§) becomes 



A tanh X = h 
5 



:i2) 



and we see that for positive h there is one bound state [jT^, A = Af,. When h is close to 
zero its value expands in powers of h 

3 45 



and hence one has 



i-l + 5 + «"" 



It is important to note that for positive h expressions (|T0|) and (|TT]) obtained by using 
the Mittag-Leffler theorem include the contribution of the bound state |T^. The small 
h behaviour of (|T0|) and (0) is due to the bound state as one can see by comparing 
(p!0|), (pH]) and (|13]). In particular, this explains why the expression ( [T0| ) is negative when 
< /i<2. 

The C-function we want to compute is defined for the part of the spectrum with uj"^ > 0. 
Hence we have to exclude the bound state. Subtracting ^ from ( ]TT| ) we get 

i?-4C(4) = y ^ = — + —h - —h^ + O(h^) . (14) 

^^ ' i^^a^ 80 945 675 ^ ^ ^ ' 



The divergence of the sum pT| ) at small h is due to the bound state. After the subtraction 
the sum becomes finite aX h = +0. In fact, it is bounded for all /i > 0, monotonically 
increasing from -^ at small /i to | for infinitely large h. We conclude that for positive h 
the function /?~^^(4) < | and the Bekenstein bound perfectly holds. 

Consider now the case of negative h. There is no bound state in this case, so that the 
expression (|TT|) gives us exactly -R~^C(4). As function of h it approaches | at /i ^ — oo 
and grows as (^ — J^) when h is close to —0. Thus, there is no upper bound for C(4) and 
the universal entropy bound can not hold in this case. 

It is easy to understand why this happens. The lowest root of the equation (0), which 
for ho = 0, hi = h reads 

atana = —h , (15) 

is always between and 7r/2 when h < 0. For small negative h 

al = -h- hi-" - ^e + 0{h') (16) 



it approaches zero and when h = it becomes the known zero-mode of the Neumann 
boundary value problem. A higher eigenvalue a„, n > 1 lies between ^{n — 1) and |n 
for all values of h. The value of C(4) for small h is then mostly due to cti. Indeed, 
approximating 

we find agreement with (^) up to 0{h^) terms. The C(4) is unbounded (that indicates a 
violation of the entropy-to-energy bound) because the lowest energy level oji = aiR~^ of 
the spectrum can be made arbitrarily small. 

For positive h there is no state with energy close to zero, the lowest a appearing in 
the interval between 7r/2 and vr. When h passes from the negative to positive values the 
lowest excited one-particle state (with eigenfrequency Ui = aiR~^) becomes the bound 
state (with u^ = —XlR~^) and should be excluded. For positive h, this saves the entropy 
bound. However, the possibility to make the lowest energy level arbitrarily small is fatal 
for the validity of the bound when h is negative. 



3 The 3D ball 

In flat (3+l)-dimensional space-time consider a massless scalar field confined to a spherical 
cavity of radius R. The boundary condition in this case is 

( d h \ , ^ 

where r is the radial coordinate. The condition (P!7| ) has the form (^ with ^ = —jh 
since for the sphere the extrinsic curvature is i^ = -|.The Dirichlet boundary condition 
corresponds to infinite h. The energy eigenfunction 0^^ = f^^{r)Yi^n{d,ip)e^'^^ expands in 
terms of the spherical harmonics Yj^„(^, y?), the degeneracy being (2/ + 1). The equation 
on the radial function reads 

Solutions to this equation should satisfy the boundary condition ( P^Tj ) and be regular at 
r = 0. There are three types of such solutions. 



1. Zero-modes 

When a; = two solutions are possible, /i(r) = r' and f2{f) = r~'~^, only the first 
being regular at r = 0. The function /i(r) satisfies the boundary condition (|1^ only if 
the parameter h in (|1^ is some (non- negative) integer, h = Iq. Then the zero mode exists 
for / = Iq and fi{r) = r'". For h = this solution is the known zero- mode of the Neumann 
problem. The similar zero-modes appear every time when h is a positive integer Iq. In 
total, there are (2/o + 1) of them. 

2. Propagating modes 

For tu^ > the solution regular at r = is 

Ur)=N^r-'/'j,^.{ujr) , 

where A^^ is normalisation constant. Since this solution should satisfy the boundary 
condition ( ]T7| ) the spectrum is discrete u!n,i = an,iR~^, where {on,;} are the roots of the 
equation 

J;, 3(0;) 

h = l-a^^ . (19) 

We have learned from the (1 + l)-dimensional example considered above that we have to 
watch for the energy level which may be arbitrary close to zero. For a fixed / the point 
a = is the point where the function staying at the right hand side of eq.(|l9D takes the 
maximal value equal /. For negative h the lowest root which may be close to zero is ai^o 
corresponding to / = 0. One has ai^i > 2 for / > 1. When h becomes positive but less 
than 1 such root appears at / = 1, ai^i. In general, for positive h lying in the interval 
lo — l<h<lQ, where /q is positive integer, the lowest root is Oi^/q corresponding to / = Iq. 
For small {Iq — h) > one finds 

3. Bound states 

These are the regular at r = solutions with a;^ < 0, 

/A(r) = iV,r-V2/,+ .(Ar) , (21) 



where A should be determined from equation 






The function of A staying at the right hand side of eq.(p2D takes at A = its minimal 
value equal / and grows monotonically as linear function for large A. The solution to the 
equation (^2]), thus, exists only when h is positive and for a given / there may be no more 
than one such solution. If h lies in the interval lo ^ h < Iq + 1, where Iq is non-negative 
integer, solution to the equation (^) exists for / = 0, ..., Iq. The total number of the bound 
states (taking into account the degeneracy due to angles) is J2iLoi'^l + 1) = (^o + 1)^- For 
small {h — Iq) > one has 



|,,._,„,,|±|,._,,,_J|±J_, 



Xl=2ilo + ^)ih-lQ)+'-^^{h-lor+ ^^ IZn, ,A h-kf + 0{h-kf . (23) 



,„-.,.o^2^V'^ ^"^^(/o + l)^ " (/o + i)^( 
We see from eqs.(pUD and ( [23| ) that when h passes through the point h = Iq the propagating 

:'0 



state with the lowest eigenfrequency uf = R ^af , becomes a bound state with u^ 



Evaluating sums of inverse powers of the roots we can again employ the Mittag-Leffler 
theorem. At fixed / one has 

___ ^2 



'^'+^^i^.+iW= 7.11(1--) ' (24) 

a>0 " 

where 7^ is some constant (see |jl5|) and {a} are the roots (|19|), Fi^i{a) = 0, where 

Fi^iz) = zJ^iz) + (/ - h)J,^,iz) . (25) 

Up to terms of order ^'+2+^ it expands as follows 

'+*' ' 2'*iT(l + i)\ m + l) (l-h) ^ 32(! + §)(( + 1) (l-h) "j'*J 
Taking the logarithm of (^), expanding in powers of z"^ and using (^) one gets 
V ^ = 1 ( ^ {l-h + 2Y _ 1 {l-h + A) \ 

h<i i6U/ + r (i~hy (/ + §)(/ + 1) {l~h) ) ^ > 

for the sum at fixed /. In order to evaluate the (^-function we have to sum over all possible 

/ and take into account the degeneracy 

00 -| 

^-t(4) = E(2^ + l)E -r- ■ (28) 

1=0 " n,l 



We remind that the (^-function is defined for the spectrum with a^ > 0. The zero-modes 
and bound states are, thus, to be excluded in (ESl). 



When h is negative there are no bound states or zero-modes. Substituting (|27|) into 
(p^) the sum over / can be computed exphcitly and is derived in terms of psi-function as 
follows 

^^ ' 3 (2/1 + 3)3 (2/1 + 3)2 ^ ' ' 

16 . . ,s 8 , , 7r2 (2/i - 1)2 

-(^(-/i)-- + 7 + 21n2)-— )— — -(- . (29) 



(2/i + 5)(2/i + 3)3' ' ' 3 ' ' 16(2/1 + 3)2 

It seems that the function ([29|) has a pole at /i = —3/2 and /i = —5/2. However, it is easy 
to check that it is regular at those points. The function (pOf ) monotonically increases from 
the value ^ — j^a.th = —oo (Dirichlet boundary condition) to infinity when h approaches 
zero. We conclude that the zeta-function is unbounded near h = 0. The expression ( [29| ) 
is valid also for the positive h. However, in this case there may be bound states and 
zero-modes (their appearance signals in that ( [29| ) has poles at integer h) the contribution 
of which must be subtracted from the right hand side of (p9|). 



When h approaches from below any non-negative integer /q the expression (|2^) for 



/ = /o grows to infinity. It is due to the fact that the lowest root (PPD, «!,/(,, becomes 



arbitrarily small. Indeed, approximating the sum ( P7| ) for / = Zq by l/<yf /^ we find that 

OO 1 1 1 

n=l "n,io "l,«o ^^y'-O + 2) ^'0 + 2) 

Thus, for h less but very close to /q one finds that the ^-function 

R-'m = ^^^^+oiio-hr , (31) 

where a^i^ is given by eq.(^), is not bounded from above. This also can be seen from 
the analysis of the poles in (^). 

The whole picture changes when h passes through the point h = Iq and becomes 
slightly grater than Iq. In this case the root ai^i^ disappears. But, instead, there appears 
a bound state, A2^. In this case, as was explained in the previous section, the expression 
(p^) for / = Iq (and, hence, also the right hand side of eq.(^)) contains the contribution 
of the bound state which should be excluded |]T9[. The sum X^n l/c^n/o ^^en is over the 



10 



propagating modes only and is finite when h approaches /q from above. As a result, the 
^-function (pHf ) is finite when h —>■ Iq + 0. A special case is when h = Iq exactly. In this 
case there appears a zero-mode /i = r^° which should be excluded when one computes the 
(■-function. When h increases, this zero-mode becomes the bound state corresponding to 
/ = /q while the rest of the spectrum changes continuously. Therefore, the case h = Iq can 
be achieved by taking the limit h —* Iq + 0. All this repeats every time when h becomes 
close to some positive integer. The general behaviour of the spectrum and the (-function 
near every such point is similar to that we had in (1 + l)-dimensional case. We conclude 
that in a way similar to the two-dimensional case the four- dimensional ((4) considered 
as function of the parameter h in the boundary condition ( p!7D is not bounded near each 
point h = Iq at which the radial equation (ITB]) has a zero-mode. 



4 The massive field 

In this section we briefly discuss the massive field. In the presence of mass m = fiR~^ the 
eigenfrequencies Ua are defined as R^uj"^ = jji^ + a^, where {a} are the roots considered in 
section 2. The (-function then reads 

^"'C(4)= E T-Ar-T^-2 ' (32) 

where sum is over states with a;^ > 0. We see that due to mass some of the bound states 
of the massless field become now the propagating states. In particular, the zero-mode is 
the state (|2T|) with X = fi. The partial summation in (|3^) for fixed / can again be done 



with the help of the Mittag-LefHer theorem |jT5|]. After analytical continuation 2; — > «/i 



eq.(p^ reads 



/.-('+^)$,^.(/.) =7.11(1 + ^) ^ 

$,+1 (/i) = /i/;^ 1 (/i) + (/ - h)ii^, (/.) . (33) 



Taking the logarithm of (^) and differentiating with respect to fi one finds 

11 



Of course, using (|^) for evaluation the sum (^) we should watch for the bound states 
(/i^ + a^ < 0) and zero-modes (/x^ + a^ = 0) and subtract their contribution from the 
right hand side of (|3^). We have two parameters in our disposal, h and /i. When h is 
negative there are no states with /x^ + a^ < 0. Let us fix /i, l^ < h < Iq + 1, where /q is 
some non- negative integer. In this case, there are (/q + 1) roots A^^ < ... < Ai < Aq of the 

A zero-mode of the massive field appears when /i equals to one of these A, say /i = A^. 
For 11 close to A^ one finds that 

where Ch is some constant. 



Substituting this expression into (|3^) we find that the sum 



diverges when /i approaches A^. If /i is slightly grater than A^, the A^-state is a propagating 
state. In fact it is the state with the lowest positive eigenfrequency, uj = fi'^ ~ A|,. The 
expression (^) is divergent because the lowest energy level ui can be arbitrarily small. 
This results in the corresponding divergence of the zeta-function (p2D. 

Decrease now /i so that it becomes slightly less than A^. Then the A^-state becomes a 
bound state of the massive field. Subtracting its contribution from ( P5| ) and ( P^ ) we get 



a finite expression. This means that the C- function is bounded from above as /i — > A^ — 0. 
This resembles the picture we had in the previous section when parameter h approached 
a positive integer. In both cases the C-function is unbounded in half-vicinity of a point 
in the parameters space where a zero-mode appears. Also, it is so because the lowest 
eigenfrequency ui is arbitrarily small when that point is approached. 
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5 The estimate for "forbidden region" 

It follows from the analysis just given that the bound (|lD can not hold for all values of 
the parameter h in the boundary condition (|17|). However, since ^(4) as function of h is 
sharply picked at integer values of /i, one may expect that the bound (|I|) holds for almost 
all values of h except a narrow region near integer h. Then (^ would be still valid for 
a "typical" boundary condition. It is indeed the case. In order to show this we need to 
estimate the "forbidden region" in the space of boundary conditions where the bound is 
violated. In order to make the estimate as precise as possible it is the best to use the 
formula (^ and its generalisation |TB[ for arbitrary degeneracy. Using (^) we find for h 



close to an integer n that 

uy^^{S{E)/E)^R ^^^- , (36) 

^(2n + 3)(n-/i) 

where ^f^ = 2ri + 2. The bound (|l|) is exceeded for h lying in the interval n — en ^ h < n, 
where for e„ we find using ( |36D the following estimate 

1 i^ngnf 



(37) 



"" 47r2(2n + 3) ' 

In particular, we find that cq ^ 0.00406, eio ^ 0.0105 and eioo — 0.00352. We see that the 
"forbidden region" near each integer n is indeed very narrow and shrinks for large n. It 
follows that the bound (|1|) holds for almost all negative values of the parameter h except 
the narrow region — eo < h < 0. Recalling the relation between h and the non-minimal 
coupling ^ in (^, /i = —4^, we find that the forbidden region for positive ^ is from to 
10^'^. Of course, the conformal coupling ^ = | is quite far from this region. 

For positive h there is a piece of the forbidden region near each integer value of h. In 
order to estimate how dense this region is in the parameter space, let us consider large 
interval < h < N. The relevant density then is given by quantity j^ J2n=i ^n- Using (p7| ) 
we find 

^|:^..-^^(>.>A')' . (38) 

for large A^. So, picking randomly a boundary condition with a positive h we almost 
always put finger on the right one for which the bound (|I|) comes ok. For example, for 

13 



the interval from to A^ = 100 the probabihty to choose a wrong condition is 4 x 10 ^. 
Thus, the bound (|I]) still holds for a "typical" boundary condition. 

6 Concluding remarks 

Recall the basic assumptions (see P): i) S{E) is microcanonical entropy defined as log- 
arithm of the number of accessible states with energy up to E; ii) E is the energy over 
vacuum (zero- modes are excluded); iii) all interactions are negligible except the boundary 
interaction expressed as a boundary condition; iv) walls of the cavity are ignored. Among 
these assumptions, the last one is perhaps the most suspicious although we just follow the 
prescriptions of paper |^. By the original idea of Bekenstein (recently re-stated in |TIi|) 
the bound (^ applies to a complete system with E being the total gravitating energy of 
the system. So it might be that when the boundary condition of the form (^ is imposed 
the boundary itself may carry some part of the total energy and, possibly, entropy. The 
bound d^) then should apply to the complete system of the quantum field and walls. This 
possibility should be further investigated. 

Alternatively, one may try to redefine the ground state. Consider for simplicity the 
two-dimensional field system analysed in section 2. The lowest energy level can be defined 
as a new vacuum, then the next excited level has energy {002 — cui) which is bounded from 
below for all values of h. So that the entropy bound holds in this case. However, the 
gravitating energy is not something one can define by a random choice of zero. It is not 
clear at the moment how the gravitating energy is actually computed for the system under 
consideration. Although, there still might be possible, by changing the rules appropriately, 
to save the bound for all boundary conditions we prefer to interpret our results, provided 
the bound applies to the field only, as imposing certain (quite relaxed) constraints on the 
value of the boundary interactions. 
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is A^ = —hohi and the roots of eq. (^ are a„ = 1112,12 = 1,2,.... Including both the 
bound state and the states with a^ > it is straightforward to evaluate the sum 
— ^ + J2'^=i ^ = h\~ + 6' ■'-^ exactly coincides with what we get by expanding the 
logarithm of eq.(|^) in powers of z'^. The contribution of the bound state is evident. I 
thank J.S. Dowker for this example. 

[19] One should also subtract the contribution of all other bound states appearing for 
/ = 0,..,/n-l. 
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